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Abstract. I review recent progress in relating meson spectral function to 
imaginary time correlation function at finite temperature calculated on isotropic 
as well as on anisotropic lattices. Special attention is payed for the lattice artifacts 
present in calculation of meson spectral functions. Results in the case of light 
quarks as well as heavy quarks are reviewed which indicate in particular that 
even in the chiral limit meson spectral functions have non-trivial structure and 
the ground state quarkonia survive up to temperature 1.5T C . 
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1. Introduction 

Spectral functions of mesonic operators play an important role in finite temperature 
QCD. Many experimental results in high energy heavy ion collisions ( e.g. low invariant 
mass dilepton enhancement, anomalous J/tp suppression etc.) can be understood in 
terms of medium modifications of meson spectral functions [[] . 

It is commonly believed that lattice QCD is capable only for calculation of static 
quantities at finite temperatures, such as the transition temperature, equation of state, 
screening lengths etc. However, it was shown by Asakawa, Hatsuda and Nakahara 
that using the Maximum Entropy Method one can in principle reconstruct also meson 
spectral functions. The method was successfully applied at zero temperature [30] and 
later also at finite temperature El El HD OH ID] ■ Though systematic uncertainties 
in the spectral function calculated on lattice are not yet completely understood, it 
was shown in Ref. [5] that precise determination of the imaginary time correlator can 
alone provide stringent constraints on the spectral function at finite temperature. 

In this contribution I am going to review recent results on meson correlators and 
spectral functions in finite temperature QCD. The rest of the paper is organized as 
follows. In section 2 I will discuss the relation between the imaginary time correlators 
and meson spectral functions. In section 3 numerical results for the light quarks will 
be discussed. Section 4 contains results on charmonia at finite temperature. Finally 
conclusions are given in section 5. 
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2. Meson correlators and spectral function in continuum and lattice QCD 

In finite temperature field theory one usually considers two types of real time 
correlation function of some operator O 

D>(t,t') = (d(t)6(t')) T , D<(t,t') = (d(t')d(t)) T , (i) 

where (...)t = {■■■e~ B '^ t )t denotes the thermal average The imaginary time 

correlation function which can be calculated also by using lattice simulations is defined 
as simple analytic continuation 

G(r) = (TO(-ir)6(0)). (2) 

( T stands for time ordered product). The spectral function can be defined through 
the Fourier transform of D > ( < \t) as 

, i D>(uj) ~ D<(uj) 1 . . 

a w) = ^- ^ = -ImD R (cu), (3) 

Z7T 7T 

where Dr(uj) is the retarded correlator. Using Eq. (J2J and KMS condition on D > ( < > 
|12| one can easily derive the following integral relation between the imaginary time 
correlator and the spectral function 

„r s f°° , , x cosh(cj(r - 1/(2T)) f°° i \ / \ / \ 

G(t)= / duja(uj) . , J -= / dua(uj)K(uj,T) (4) 

Jo smh Jo 

Using this relation one can in principle reconstruct the spectral function by calculating 
G(t) on lattice. At zero temperature the kernel, K{oj, r) reduces to simple exponential 
and at large Euclidean times the correlation function picks up the contribution from 
the lowest lying meson state in o~(co), i.e G(r) = exp(— mr). At finite temperature 
the analysis of the asymptotic behavior of the correlators is no longer possible as r 
is limited to the interval [0, 1/T] where excited states are equally important as the 
ground state. Additional complications arise in lattice calculations where correlators 
are calculated only at finite set of Euclidean times tT = k/N T , k = 0, ...,N T — 1 
with N T being the temporal extent of the lattice. In order to reconstruct the spectral 
functions from this limited information it is necessary to include in the statistical 
analysis of the numerical results also prior information on the structure of cr(u>) (e.g. 
such as <t(u>) > for u> > 0). This can be done through the application of the 
Maximum Entropy Method (MEM) f. 

In order to study meson properties at finite temperature appropriate choice of 
the operator O should be made. One possible choice is local meson operator bilinear 
in quark-antiquark fields (current) EJ El IE] 

Oh(t, x) = Z H q(T,x)T H q(T, x), ^ = 1,75,7^,757^ (5) 

for scalar, pseudoscalar, vector and axial vector channels correspondingly. The 
normalization constant Zh relates the current calculated in lattice regularization 
scheme to current in MS scheme with ^jjg — a -1 (with a being the lattice spacing). 
We can define then the temporal correlators at finite spatial momentum p 

G H (r,p) - (0 H (r,p)O f H (r, -f)}, 0{r,p) - ]T e^O H {r,x) (6) 

X 

X Other methods of introducing prior information into the statistical analysis have been also discussed 
in Ref. I A for the zero temperature case and in Ref. ft 41 for finite temperature QCD. 
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One can also consider correlators of extended operators defined as ^3 

6 H (t, x) = ^2 0(y)9( r > x)Tq(T, x + y). (7) 
y 

Here 

m = zM-bW) (8) 

is the trial wave function which controls the size of the meson source and can be 
regarded as a physical input to the problem. The use of extended operators makes the 
reconstruction of the spectral function easier as the correlator is dominated by a single 
peak in the spectral function |10j . However, the corresponding spectral function is not 
related to a physically observable quantity, it can provide information only about 
the mass and the width of the resonance but not about the corresponding decay 
constant. This also only works for sharp resonances. At high temperature "mesons" 
will appear as broad structures in the spectral function and the above method is no 
longer applicable. 

Though the in-medium properties of mesons are encoded primarily in the 
temporal correlator and the corresponding spectral function, in finite temperature 
QCD it is customary to study spatial correlators defined as 

G H = (O(*)O+(0)>, 0{z) = Y, 0(r,x,y,z). (9) 

Since the lattice extent is not limited in the spatial directions one can study their 
large distance behavior. At large distances the spatial correlators decay exponentially 
and the exponential decay governed by the so-called screening mass. The spatial 
correlators can be related to the spectral function via following relation 

d Pz e^ z G(t,Q,0, Pz )= dp,e*«« / ^ 1 ' ' ' PzJ (10) 

-oo JO J-oo JO U 

In general the screening masses are different from the pole masses, however, from the 
above equation one can easily see that in the special case when the spectral function 
is dominated by a single ^-function for small u) the screening and the pole masses are 
equal. 

The relation between the imaginary time correlator and the spectral function (0} 
is valid only in the continuum theory. As the correlators are calculated on lattice the 
question arises whether the spectral representation of the form (0} can be derived also 
for the lattice correlator. In ref. JS] it was shown that this is indeed the case in the 
free theory and that the cutoff effects present in the correlator at small separations 
are contained in the spectral function. Because of the asymptotic freedom calculation 
in the free theory are relevant at high temperature and/or large energies ui. In Fig. 
^1 show the ratio of of the lattice spectral function to the corresponding continuum 
spectral function versus au; as well as the lattice pseudoscalar spectral function versus 
uj/T in the free theory. As one can see from the figures the lattice spectral functions 
start to deviate from the continuum one for Lua > 0.5 and in the pseudoscalar and 
vector channels show a peak like structure around uia ~ 2. The lattice spectral 
functions vanishes for loo > ui max a ~ AN T . As the temporal extent of the lattice N T is 
increased for fixed temperature, i.e. the lattice spacing decreases (a — 1/(N T T)) the 
peak structure moves to larger value of u>/T (see Fig. [Q. Peak structure around uia ~ 2 
was also observed in the interacting theory 0]. The main reason for the large cutoff 
effects in the lattice spectral function is that on lattice the quark dispersion relation 
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Figure 1. The ratio of the lattice and continuum spectral functions in different 
channels (left) versus uja and the lattice pseudoscalar spectral function versus 
lu = u)/T for different values of N T calculated using Wilson fermions. See Ref. 
16 for further details. 



gets modified at large momenta, pa > 1, relative to its continuum counterpart. Cutoff 
effects in the spectral function can be substantially reduced by using an improved 
fcrmion action, the so-called truncated perfect action on hypercube |T7 for which the 
lattice quark dispersion relation is much closer to the continuum even for pa > 1 |16| . 

3. Numerical results in the light quark sector 

Meson correlators and spectral functions for small quark masses were intensively 
studied during recent years [SJ EH ■ In Kefs. |SJ E] spectral functions were 
studied using non-perturbatively 0(a) improved Wilson action J5| and isotropic 
32 3 x 16, 48 3 x 12, 64 3 x 16 and 64 3 x 24 lattices. In these studies quark masses 
corresponding to the pion masses in the range 400MeV to lGev were considered in the 
confined phase, while in the deconfined phase calculation were performed in the chiral 
limit and the temperature interval was 0AT c — 3T C with T c being the deconfinement 
temperature. The renormalization factors Zh appearing in Eq. were determined 
non-perturbatively in Ref. |19| for the vector and axial vector channels, while for 
the scalar and pseudoscalar channels they were calculated using tadpole improved 
perturbation theory. Asakawa, Hatsuda and Nakahara have studied meson spectral 
functions using 32 3 x iV r anisotropic lattices with a T /a a = 4 and N T = 96 — 32 
corresponding to temperatures 0.8T C — 3T C [§]■ The quark masses used in this study 
correspond to m w /m p = 0.7. Below the deconfinement temperature neither the 
correlators nor the spectral functions indicate in-medium change of meson properties 
contrary to existing theoretical predictions pQ. This is likely due to the quenched 
approximation and large quark masses used in these studies. 

Before discussing the results on the spectral function above deconfinement let 
us first look at the correlators. The correlation functions alone can provide some 
constraints on the spectral functions, for example the analysis of the vector correlation 
function provides stringent constraint on the thermal dilcpton rate At high 

temperature one would expect that the quark propagators should be close to the 
free ones. Therefore in what follows we will always normalize the meson correlators 
by the corresponding lattice correlators in the free theory. We will restrict the 
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Figure 2. The ratio of meson correlators at 1.5T C in different channels to the 
corresponding correlators in the free theory calculated on 64 3 X 16 lattice. 
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Figure 3. The ratio of meson correlators at 1.5T C and 3T C in the pseudoscalar 
(left) and vector (right) channels to the corresponding correlators in the free 
theory calculated on lattice. 



discussion to the case of zero spatial momentum p = 0. In Fig. [21 I show the 
meson correlators at 1.5T C in different channels calculated on 64 3 x 16 lattice. As 
one can see from the figure the scalar-pseudoscalar and vector-axial- vector correlators 
become degenerate except at very small tT where lattice artifacts present in the Wilson 
formulation explicitly break chiral symmetry. This can be viewed as an indication of 
the Ua{^) symmetry restoration at high temperature. In Fig. |21l show the vector and 
pseudoscalar correlators for 1.5T C and 3T C for different values of N T . From the figure 
it is evident that the vector correlator stays close to its free value even at 1.5T C while 
the pseudoscalar correlator is considerably enhanced compared to the corresponding 
free correlators. Note also that in the vector channel the correlators is N T (i.e. lattice 
spacing) independent, while in the pseudoscalar case there is a small N T dependence 
which is probably due to the perturbative error in calculation of Zps- In Fig. 01 1 
show the spectral functions for the pseudoscalar and vector channels. The spectral 
functions in both channels show a peak like structure at lu ~ (5 — 6)T which is more 
pronounced in the pseudoscalar case and probably leads to the large enhancement of 
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Figure 4. The spectral functions in the deconfincd phase for the pseudoscalar 
channel (left) and the vector channel (right) reconstructed using MEM on 64 3 X 16 
lattice. 




5 10 15 20 25 30 5 10 15 20 25 30 

ra[GeV] <o[GeV] 



Figure 5. The spectral functions p{u>) = cr(ui)/ui 2 in different channels calculated 
at 1.4T C (left) and at 1.9T C (right) using anisotropic lattices [3] Note that 
normalization of the spectral functions is different from the case on isotropic 
lattices. 



the pseudoscalar correlator over the free case discussed above. Note that the position of 
the peak appears to be proportional to the temperature. Spectral function calculated 
on anisotropic lattice |3] are shown in Fig. [5] and exhibit a similar peak like structure, 
roughly for the same value of oj/T. However, at lower temperature, T — 1AT C the 
corresponding structure in the spectral functions appears to be more sharp (see Fig. 

Meson correlators in the spatial directions and the corresponding screening masses 
have been studied since long time both in full and quenched QCD [2U1 [2T1 WR I2TH 1251 
|^] [2EI 123 HH • However, reliable determination of the screening masses is available 
only in quenched approximation |241 1271 12"5] as lattices with temporal extent N T > 8 
are needed in such analysis |27| . The most recent results for the screening masses are 
summarized in Fig. |SJ The screening masses are close to their asymptotic (T — > oo) 
value 2ttT already at temperature for T = 1.5T C . This could lead to the conclusion that 
there is an almost free quark propagation in spatial directions in the deconfined phase 
at temperatures as low as 1.5T C . However, a closer look on the spatial correlators in 
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Figure 6. The screening masses (left) and spatial correlators (right) for vector 
and pscudoscalar channels at T > T c . The screening masses have been calculated 
in quenched QCD using improved Wilson fermions 1241 1281 and Kogut-Susskid 
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Figure 7. The ratio of spatial and temporal correlators in the deconfined phase 
to the corresponding free correlators for the pseudoscalar (left) and vector (right) 
channels. 



Fig. EJreveals that this is not the case even at 3T C . The correlators calculated by using 
lattice simulations are several times larger than the corresponding free ones. This is 
not unexpected since it is well known that the physics at distances z > 1/ (g 2 T) is 
non-perturbative even at very high temperature |29| . In Fig. 01 show the vector and 
pseudoscalar spatial correlators at short distances together with the temporal one. 
As one can see from the figures at small separations temporal and spatial correlators 
show very similar behavior. 

4. Heavy quarkonia correlators and spectral functions 

Heavy quarkonia correlators and spectral functions were studied in Refs. [101 I15| 
using anisotropic lattices and extended operators as well as on isotropic lattices and 
point correlators In Ref. |S] pseudoscalar, vector, scalar and axial vector channels 
were considered which correspond to 1 Sq (t] c ), 3 <Si (J/"4>), 3 Po (Xco) and 3 Pi (xd) 
charmonia states respectively. As in the light quark sector some statements about 
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Figure 9. Spectral function in the pseudoscalar (left) and axial vector channels 
(right) reconstructed using MEM. 



in-medium properties of charmonia can be made just by analyzing the behavior of 
the correlators. In Fig. 01 show the ratio of different quarkonia correlators to the 
corresponding free correlators. As one can see from the figures the correlators in the 
pseudoscalar and vector channels show very little change across T c while in the case of 
the scalar and axial vector channels larger changes in the correlators are visible. This 
implies that the ground state charmonia ( 1 5'o and 3 S\) are very likely to survive in the 
deconfincd phase. The spectral function reconstructed using MEM are shown in Fig. 
Elwhich shows indeed the the ground state peak survives in the deconfined phase while 
the excited 3 Po state is dissociated already at T — 1.25T C . Spectral functions were also 
reconstructed for 1.5T C and and show again that the ground state charmonia survive 
while the excited P-states are dissociated. Ground state charmonia ( 1 S'o, 3 Si states) 
spectral functions were also studied in Ref. using extended operators. The results 
for the vector channel are shown in Fig. The J j ip mass (the peak position) changes 
very little across T c while the corresponding peak in the spectral function gets broader 
above T c . This broadening of the ground state peak is also visible in calculations with 
local operators. However, such broadening is not necessarily a physical effect and may 
come from the fact that less data points N T are available at higher temperature [3] 
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input wave functions (different values of b). 



and also the temporal extent of the lattice in physical units become smaller ■ The 
authors of Rcf. have found that the width of the peak in the spectral function also 
depends on the trial wave function <fi(y) used in constructing the extended operator 
(see Fig. I10|) . the peak is broader for larger value of b in the trial wave function 4>(y). 
Therefore MEM cannot yet provide a reliable estimate of the width of the J/ -0 in the 
plasma. To overcome this difficulty \ 2 -fits to the Breit-Wigner form of the spectral 
function was used |10|. This led to a first hints for a non-zero widths of charmonium 
states above T c . 

The existence of ground state charmonia at 1.5T C is in sharp disagreement with 
predictions of potential models based on color screening 301 EH A possible reason for 
this discrepancy is an oversimplified picture of color screening used in these studies. In 
fact recent analysis of the free energy of quark-antiquark pair suggest quite complicated 
medium dependence of inter-quark forces |H2I 02] . 

5. Conclusions 

In this paper recent results on meson spectral functions and correlators calculated on 
lattice were discussed. Though systematic uncertainties are present in the spectral 
functions the results appear to be very interesting and even intriguing to some extent. 
The presence of the peak like structures in the meson spectral functions well above 
deconfinement temperature is, too large extent, unexpected. These peak structures 
seems to be present in the spectral functions calculated on isotropic lattices [5] |S] 
as well as on anisotropic lattices [SJ and is some case their presence is supported 
by direct inspection of the corresponding correlation functions (c.f. Fig. [3] (left) 
and Fig. |HJ. In the case of heavy quarkonia the interpretation of the peaks in the 
spectral function is clear, they correspond to the quarkonia states which survive in the 
deconfined phase up to temperatures 1.5T C . We have seen that the position of these 
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peaks does not change considerably compared to the zero temperature case. The 
interpretation of peak-like structures in the light quark sector is less evident as the 
position of these peaks seems to be proportional to the temperature. Certainly much 
more work is needed in order to be able to interpret the present findings in physical 
terms as well as for the detailed understanding of the systematic uncertainties involved 
in extraction of the spectral functions from lattice correlators. 
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